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Abstract 



Using more than 30 years of lunar laser ranging (LLR) data, the lunar orbit has been modeled 
and fit with several millimeters precision. This fit provides comprehensive confirmation of the 
general relativistic equations of motion for solar system dynamics, achieving several key tests 
of Einstein's tensor theory of gravity, and strongly constraining presence of any additional 
long range interactions between bodies. Earth and Moon are found to fall toward the Sun at 
rates equal to better than a couple parts in 10 13 , confirming both the universal coupling of 
gravity to matter's stress-energy tensor, and gravity's specific non-linear coupling to itself. 
The expected deSitter precession (with respect to the distant 'fixed' stars) of the local inertial 
frame moving with the Earth-Moon system through the Sun's gravity is confirmed to 3.5 
parts in 10 3 precision (± .07 mas/year), and Newton's gravitational "constant" indeed 
shows no cosmological time variation at the part in 10 12 per year level. Most all of the 
1/c 2 order, post-Newtonian terms in the N-body equations of motion — motional, gravito- 
magnetic, non-linear, inductive, etc. — contribute to the measured details of the lunar orbit, 
so LLR achieves near-completeness as a gravity experiment and probe. The precision of these 
measurements, especially those connected with lunar orbit frequencies and rates of change 
of frequencies, will further improve as LLR observations continue into the future with use of 
latest technologies. 
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1. Introduction. 

The precise fit of the lunar laser ranging (LLR) data to theory yields a number of the most 
exacting tests of Einstein's field theory of gravity, General Relativity, because almost any 
alternative theory of gravity predicts a number of changes (from that produced by General 
Relativity) in the lunar orbit which would be readily detected in the LLR data. Some 
of the most interesting and fundamental of such theory- dependent effects, and which are 
particularly well-measured by LLR, include 1) a difference in the free fall rate of Earth and 
Moon toward the Sun due to gravity theory's non-linear structure acting on the gravitational 
binding energy within the Earth, 2) a time variation of Newton's gravitational coupling 
parameter, G — > G(t), related to the expansion rate of the universe, and 3) precession of the 
local inertial frame (relative to distant inertial frames) because of the Earth-Moon system's 
motion through the Sun's gravity. 

Measurements of the round trip travel times of laser pulses between Earth stations 
and sites on the lunar surface have been made on a frequent basis ever since the Apollo 11 
astronauts placed the first passive laser reflector on the Moon in 1969. Today about 15,000 
such range measurements are archived and available for use by analysis groups wishing to 
fit the data to theoretical models for the general relativistic gravitational dynamics of the 
relevant bodies, the speed of light function in the solar system, tidal distortions of Earth and 
Moon, atmospheric corrections to light propagation, etc. An individual range measurement 
today has precision of about a centimeter (one-way), but a new generation observing program 
plans to improve this range measurement precision down to a millimeter. Because of the 
large number of range measurements, some of the key length parameters which describe the 
lunar orbit are already estimated with precision of a few millimeters, and key lunar motion 
frequencies to fractional precisions of a few parts in 10 12 . 

Because both the Earth's mass and that of the Moon are sufficiently large, the orbits 
of these bodies can be modeled as single orbital "arcs" extending over three decades through 
time. The complete model used to fit the many range measurements contains in excess of a 
hundred parameters P m which are optimally adjusted from their nominal model values P$ 
by amounts 5P m = P m — P$ determined in a weighted least squares fit type procedure 

N M 

Minimize J2 W v J2 [f(m)i5P m - rj [f(n)jSP n - r,] 

i,j=l m,n=l 

with the N range measurements being identified by the labels % and j, and the M model 
parameters being identified by the labels m and n. Wij are the weightings given to each 
measurement (pair) and are usually taken to be diagonal in ij and inversely proportional 
to the square of inferred measurement errors; the residuals Ti are the differences between 
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observed and calculated range values, r\ = R b s (ti) — R C aic(ti); and the remaining functions 
f(m)i are the parameter partials which give the sensitivity of the modeled (calculated) range 
to change in each model parameter value 

r / \ dR ca i c (ti) 
f( m )i = — o75 

evaluated at the time U of the ith range measurement. 



Among the very many model parameters, the information needed for testing relativis- 
tic gravity theory is concentrated in only a handful of orbital features. The needed orbital 
parameters are connected with four key oscillatory contributions to the lunar motion, the ec- 
centric, evective, and variational motions and the parallactic inequality, which are illustrated 
in Figure 1. The eccentric motion produces an oscillatory range contribution proportional to 
cos(A), A being anomalistic (eccentric) phase, and is a natural and undriven perturbation 
of circular motion. The variation is driven by the Sun's leading order quadrupolar tidal 
field and produces a range contribution proportional to cos(2D), D being synodic phase 
from new moon. The parallactic inequality is driven by the Sun's next order octupolar tidal 
field, and its range perturbation cos(D) has monthly period. The evection is a hybrid range 
perturbation proportional to the eccentric motion as modified by the variation and having 
time dependence cos(2D — A). The eccentric and evective motions, which alter the times of 
eclipses, were discovered by the anchients; the variation and parallactic inequalities, which 
do not alter times of eclipses, were only found during and after the era of Newton. 

The amplitude of the parallactic inequality, Lpj, is unusually sensitive to any difference 
in the Sun's acceleration rate of Earth and Moon [2]. The frequency of the eccentric motion, 
the anomalistic frequency A, when compared to other lunar frequencies determines the 
precession rate of the Moon's perigee. This rapid precession, which completely rotates the 
orbit's major axis in about 8.9 years, is primarily driven by the Sun's tidal acceleration, 
but there is a leading order relativistic contribution to this precession rate interpreted as 
an actual rotation of the local inertial frame, the deSitter precession, resulting from motion 
through the Sun's field of gravity. From measurement of the time rates of change of the 
Moon's anomalistic and synodic frequencies, A and D, a rather clean measurement can be 
made of a time rate of change of Newton's coupling parameter G. The Earth-Moon range 
model can be expressed in terms of these primary contributions 

r{t) = L Q — L ecc cos(A) — L evc cos{2D — A) — L var cos{2D) — Lpicos(D) + ... 

with phases advancing as A = A Q + A(t — t Q ) + A(t — t Q ) 2 /2 + ... and similarly for synodic 
phase D. The LLR measurement of Lpj, A, D, A, and D are the foundations for the 
gravity theory tests. 
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Lunar Orbit's Four Main Perturbations 



Eccentricity (20905 Km) Evection (3699 Km) 




toward Sun 

.1 




Variation (2956 Km) 



toward Sun 




Parallactic Inequality (110 Km) 



Figure 1. Four lunar orbit perturbations from a nominal circular orbit (dot- 
ted) are shown. They produce oscillatory Earth-Moon range terms: the eccen- 
tric oscillation ~ cos(A), the variation oscillation ~ cos(2D), the parallactic 
inequality oscillation ~ cos(D), and the evective oscillation ~ cos(2D — A), 
with respective amplitudes indicated. Key tests of general relativity arc 
achieved from precise measurements of amplitudes or phase rates of these 
perturbations. Measurement of the amplitude of the parallactic inequality 
determines whether Earth and Moon fall toward the Sun at same rate. Mea- 
surements of the synodic phase D and anomalistic (eccentric) phase A rates 
and rate of change of these rates determine the deSitter precession of the 
lunar orbit and time rate of change of Newton's G. 



2. Dynamical Equations For Bodies, Light, and Clocks. 
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LLR comprehensively tests the 1/c 2 order, gravitational N-body equations of motion which 
analysis groups integrate to produce orbits for Earth, Moon, and other relevant solar system 
bodies. The Sun-Earth-Moon system dynamics is symbolically illustrated in Figure 2, with 
the rest of the solar system bodies sufficiently considered at the Newtonian level of detail. 
The Earth moves with velocity V and acceleration A with respect to the Sun, while the Moon 
is moving at velocity V + u and acceleration A + a. (If preferred frame effects were to be 
considered for cases when gravity is not locally Lorentz- invariant, the Sun's cosmic velocity 
W also becomes involved.) [4] There are a variety of post-Newtonian forces acting on Earth 
and Moon by the Sun, by each other, and on themselves (self forces) which are dependent on 
these general motions. Included in these are non-linear gravitational forces for which each 
mass element of the Earth and Moon experiences forces due to the interactive effect of the 
Sun's gravity with the other mass elements of the same body, or of the other neighboring 
body. The accelerations of individual mass elements of Earth also induce accelerations on 
the other mass elements of Earth, and similarly with the Moon. Acceleration of Earth 
induces an acceleration of the Moon. Altogether, these 1/c 2 order accelerations produce a 
rich assortment of modifications of the Earth-Moon range which LLR can measure. 

The N-body equation of motion in metric gravity has been formulated in the literature 
for the completely general case [16]. Not observing any violations of local Lorentz-invariance 
or breakdown of conservation laws in solar system gravity, I here specialize consideration to 
the fully conservative, locally Lorentz- invariant, lagrangian-based gravitational equation of 
motion (plus cosmological variation of Newton's G). For N bodies in general motion and 
configuration, and valid for a broad class of plausible metric theories of gravity, scalar-tensor 
theories in particular, the order 1/c 2 equations of motion for these N bodies take the form 
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Velocities and Accelerations of Sun, Earth, and Moon 



— * 

v 




Figure 2. When formulating the Earth- Moon dynamics in the solar system 
barycentric frame, there are post-Newtonian force terms acting between Sun, Earth 
and Moon which depend on either the velocity or acceleration vectors of both the 
Earth and Moon. Body self- accelerations also result from the inductive inertial 
forces acting between the mutually accelerating mass elements (i, j) within each of 
these bodies. The intrinsic non-linearity of gravity also produces net external forces 
on these bodies proportional to not only the presence of other bodies, but also to 
their internal gravitational binding energies. The motional, accelerative, and non- 
linear contributions to the three body system's dynamics, taken collectively, make 
LLR a comprehensive probe of the post-Newtonian dynamics of metric gravity in 
the general case. If the dynamics is not locally Lorentz invariant, then the velocity 
W of the solar system through the cosmos leads to novel forces and resulting ob- 
servable effects in LLR proportional to W (or its square); but such effects have not 
been seen. 
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F - \v% - a t ■ v t v t - (2 7 + 1) ]T ^ (1) 

with ?Tj = dfi/dt, di = dvi/dt, = \fi — fj\, and i,j,k = 1 ... N. The speed of light factor 
1/c 2 has been set equal to 1 in lines B ... F to simplify presentation. The body gravitational 
mass strengths = GM(G)i are indicated along with the Newtonian acceleration vectors 



7 and (3 (with (3* = 2f3 — 1) are two Eddington parameters which quantify deviations of 
metric gravity theory from Einstein's pure tensor theory in which both these parameters 
equal one. The several lines of this total equation of motion warrant individual descriptions 
and brief discussions. 

Line A. If the metric theory Eddington parameters 7 and (3 differ from their general 
relativistic values 7 flr = f3 gr = 1, application of the equation of motion relativistic corrections 
from lines B through F to a body's internal gravity finds that the gravitational to inertial 
mass ratio of a celestial body depends on its gravitational self-energy content [1] 

M < G > = 1 _ W - 3 - 7 ) « / 4fMg> A <Py + order 1/c' (2) 



M(I) y ^ " 2Mc 2 J \x-y\ 

Another way to view the above ratio is in terms of a spatially varying gravitational coupling 
parameter G 

G(f, t) Goo [l - (4/3 - 3 - 7) U(f, t)/c 2 

in which a body with a significant part of its mass-energy coming from its gravitational 
binding energy experiences the additional acceleration 

Odi = kt^-c VG 

oG 

with the leading gravitational energy contribution to body mass being Newtonian contribu- 
tion , „ , s 

dM _ G [ P(r)P(r') d s rd 3 r , 
dlnG 2c 2 J \f— f'\ 

When cosmological equations from a metric theory are considered, Newton's coupling pa- 
rameter G will also generally be found to vary in time in proportion to the Hubble expansion 
rate of the universe 

§~ (4/3-3-7) # (3) 



8 



The presently most precise way to measure any deviation of (3 from its general relativistic 
value is through measurement of the M(G)/M(I) ratio of Earth using LLR data. 

Line B. Gravity couples to itself, thereby producing non-linear gravitational forces 
among and between bodies. 

Line C. Just as pairs of moving charges generate magnetic forces between themselves 
in proportion to the velocities of both charges, pairs of moving masses generate gravito- 
magnetic forces between themselves. This force acts between the mutually moving Earth 
and Moon and contributes to the necessary Lorentz-contraction of the lunar orbit as viewed 
from the solar system barycenter. 

Line D. Masses in motion both produce and couple to gravitational fields differently 
than masses at rest. The package of velocity-dependent acceleration terms in this line plus 
line C lead to local Lorentz-invariance of gravity. Any further modifications of this package 
(beyond the 7-dependence) will lead to additional terms in the equation of motion with one 
or two powers of body velocities being replaced by the velocity W of the solar system relative 
to the universe preferred frame. A variety of preferred frame effects which would then result 
have been empirically ruled out in LLR and other solar system observations [5]. 

Line E. Accelerating masses generate inductive gravitational forces on other proximite 
masses. 

Line F. The inertia of a mass is altered by its motion and by its proximity to other 
masses. The combination of terms from this line plus line E are necessary in order that a 
body's gravitational self-energy contributes to its total inertial mass in accord with special 
relativity's prescription M = E/c 2 . This modification of inertia is part of the M(G)/M(I) 
calculation for a celestial body. 

LLR measures the round trip time of propagation of light between two separate body 
trajectories, and this measurement is made by a specific clock moving on a particular tra- 
jectory. So in the solar system barycentric and spatially isotropic coordinates employed to 
express the body equations of motion given by Equation (1), there are also requirements for 
the post-Newtonian modifications to the light coordinate speed function and to the clock 
rates, these respectively being 



c(f, t) 



•oo 



1 - (1+7) U(f,t)/c 2 



(4) 



and 



dr = dt 1 - v 2 /2c 2 - U(f, t)/c 2 



(5) 
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in which U (r) is the total Newtonian gravity potential function due to solar system bodies 

( 6 ) 

Because the Earth moves in the solar system barycentric frame, and it rotates at rate u, 
there must be two corrections applied to an Earth surface location a; first there is the 
Lorentz-contraction of the extended body 

8a = -a-V V/2c 2 

and because of special relativity's non-absolute nature of time simultaneity there is a further 
displacement of the rotating Earth surface locations 

5a = V ■ a {V x a)/c 2 

These light and clock equations, and special relativistic body distortion effects play only 
supportive (but necessary) roles in fitting LLR data; the main science emerges from the 
body equations of motion as given by Equation (1). 



3. LLR's Key Science-Related Range Signals. 

Assocated with each feature of gravitational theory which is tested by LLR, there are specific 
range signals in the LLR data whose measurements yield the information about theory. 
Several of these signals are here described. 

Violation of Universality of Free-fall. 

Because of gravitational self-energies (internal gravitational binding energies) in celestial 
bodies, they will generally possess gravitational to inertial mass ratios which differ from 
each other as indicated in line A of Equation (1) and given by Equation (2). But there are 
other ways in which bodies may accelerate at different rates toward other bodies. Within 
the paradigm that forces between objects are carried by a field, an additional long-range 
interaction in physical law generates a force between bodies i and j which will typically have 
the static limit form 

fi = KiVi-i-e-^ (7) 

r ij 
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The bodies' coupling strengths K t and Kj, except in special cases such as metric scalar- 
tensor gravity in which Ki ~ Mj, will be attributes of the bodies which are different than 
total mass-energy {non-metric coupling); and the dependence on distance of this force will be 
either inverse square if the field is massless, or Yukawa-like if the underlying field transmitting 
this force between bodies has mass. Such a new force will produce a difference in the Sun's 
acceleration of Earth and Moon, because the latter two bodies are of different compositions 
- the Earth has a substantial iron core while the Moon is composed of low-Z mantle-like 
materials. The fractional difference in acceleration rates of Earth and Moon amounts to 

and it will supplement any difference of accelerations resulting from the possible anomalies 
in the bodies' gravitational to inertial mass ratios due to gravitational self energies. LLR 
has become a sufficiently precise tool for measuring |5a em |, it now competes favorably with 
ground-based laboratory measurements looking for composition-dependence of free fall rates; 
and LLR is the premier probe for measuring a body's M(G) /M(I) ratio as given by Equation 
(2). 

If Earth and Moon fall toward the Sun at different rates due to either of the mechanisms 
discussed, then the lunar orbit is polarized along the solar direction. Detailed calculation of 
this polarization reveals an interesting interactive feedback mechanism which acts between 
this cos(D) polarization and the cos(2D) Newtonian solar tide perturbation of the lunar 
orbit, the variation). The result is an amplification of the synodic perturbation 

Sr(t) me = ~ RF{Vt/uj)5 me cosD (8) 

= 2.9 10 12 5 me cosD cm (9) 

with 5 em = \(a e — a m )/g s \, R is distance to the Sun, ft and u are the sidereal frequencies of 
solar and lunar motion, and D is the lunar phase measured from new moon. The feedback 
amplification factor for the lunar orbit is already F{VL/u) = 1.75; it grows further with larger 
orbits with an interesting resonance divergence for an orbit about twice the size as that of 
the Moon [14, 15]. Computer integration of the complete Eq. (1) for the Sun-Earth-Moon 
system dynamics confirms these analytically estimated polarization sensitivities. 

The most recent fits of the LLR data find no anomalies in the cosD amplitude to 
precision of 4 millimeters, so from Equation (y) 5 me is constrained to be less than 1.3 10~ 13 . 
Neglecting any possible composition dependence, and using Equation () with an estimate 
for the fractional gravitational self energy of the Earth being 4.5 10 -10 , a constraint on a 
combination of the two Eddington parameters 

|4/3-3-7| < 4 10~ 4 (10) 
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If metric gravity is a combination of scalar and tensor interaction, the small number of 
this constraint is an approximate measure of the scalar interaction strength compared to 
the dominant tensor interaction. One scenario which could explain today's weakness of the 
scalar interaction is illustrated in Figure R. Scalar-tensor metric gravity involves one coupling 
function V(<f>); the slope of this function gives the strength of the scalar interaction, and in 
combination with the function's curvature also determines gravity's 1/c 2 order non-linearity; 
near an extremum of V(<f>) the Eddington parameters are given by simple properties of the 
coupling function 

" 

As the universe expands, the dynamical equations for the background scalar field will drive 
the scalar to a minimum of the coupling function, if it exists, and where 7 and (5 take 
their general relativistic values. Scalar gravity turns itself off naturally if an "attractor" 
exists in its coupling function V ((/>). But that process, being dynamical, should not be 
entirely complete today, and the small remnant of scalar interaction may still be detectable 
by sufficiently precise testing of relativistic gravity using LLR and other experiments [8, 9]. 

The LLR result can also place limits on the spatial gradient of the fine structure con- 
stant, a = e 2 /hc, in the proximity of the Sun. If a is a function of a scalar field whose source 
includes ordinary matter, a spatial gradient of a near bodies should exist, and composition- 
dependent accelerations of other objects toward this body should occur 

dlnMi 2 Va 
olna a 

The dominant electromagnetic contribution to the mass-energy of different elements is due 
to the electrostatic energy among the Z nuclear protons. This energy fractionally varies 
by an order of magnitude (from a few parts in 10 4 to a few parts in 10 3 ) as one proceeds 
through the periodic table from low-Z to high-Z elements. For the Earth with its iron-core 
and Moon composed almost entirely of mantle-like materials, one can conclude from the LLR 
constraint on 5 me that any gradient of a due to and toward the Sun is quite small compared 
to the Sun's gravitational field g s 

° 2 igjH < 4KT 10 

\g s \ 

This should be compared with the best constraints on time variation of a, which in units of 
the Hubble expansion rate are substantially weaker 

H- < 1(T 5 

a 
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Figure 3. Typical cosmological dynamics of a background scalar field is 
shown if that field's coupling function V(<f>) has an attracting point O . The 
strength of the scalar interaction's coupling to matter, proportional to the 
derivative of the coupling function, weakens as the attracting point is ap- 
proached; so in a scalar-tensor metric theory, for example, the Eddington 
parameters 7 and (3 both approach the pure tensor gravity values of one. 
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This suggests that unless there are unusual sources for the scalar field which controls the 
value of a, e.g., sources which are present in an average cosmological context but which do 
not concentrate in ordinary matter, or other special situations, then today's LLR constraint 
on the spatial gradient of a is the significant present measure of the constancy of a. 



Geodetic Precession of the Local Inertial Frame. 



Because Earth and Moon travel at different velocities through the Sun's gravitational field, 
terms from lines D and F of Equation () are present which accelerate the Moon relative to 
Earth. A particularly interesting part of the relative acceleration is proportional to both V 

— * 

and u and the Sun's acceleration with, as shown in Figure 1, V being the velocity of the 
Earth relative to the Sun, and u the velocity of the Moon relative to Earth. These terms 
form deSitter's Coriollis-like acceleration 

5a m = 2 Q dS x u (13) 

and its geometrical interpretation is the local precession of the inertial frame at rate fids 
which amounts to about 19.2 mas/y. The effect of this perturbing acceleration on the orbit is 
primarily an additional rate of perigee precession with respect to distant inertial space. This 
is measured by comparing the Moon's anomalistic frequency A (rate of eccentric motion) 
with its synodic frequency D (rate of monthly phase), and with the latter converted into 
lunar sidereal frequency u (orbital rate) by adding to D the annual rate Q which is provided 
by results from other solar system experiments. Sidereal minus anomalistic frequency of 
lunar motion includes deSitter's precessional rate as a supplement to the Newtonian tidal 
contributions to perigee precession. These lunar frequencies are measured from range signal 
perturbations whose size grows linearly in time. The Moon's range from Earth includes 
several dominant oscillatory contributions 

5r me = L ecc cos(A) + L var cos(2D) + L evc cos(2D - A) + ... 

with L ecc being the amplitude of eccentric motion, L var being the amplitude of solar tidal per- 
turbation called variation, and L evc being the amplitude of the hybrid evection perturbation 
due to both the solar tidal force and the eccentric motion of the Moon. The least-squares-fit 
of the LLR data, which yields best estimates for the two key lunar frequencies, will then 
involve the parameter 'partials' 

d5r r 



OA 



= -t (L ecc sin(A) - L evc sin(2D — A)) 
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Q5r 

™ = -It (L var sin(2D) + L evc sin(2D - A)) 

dD 

Measurement precision of the deSitter precession grows especially with total time of the LLR 
experiment, not only because of the growing quantity and quality of the accumulated range 
measurements, but also because of the linear growth in signal sensitivity. A most recent fit 
of the LLR data confirms presence of the geodetic precession with precision of 0.07 mas/y 
[10]. 



Time Evolution of Gravity's Coupling Strength G. 

Time evolution of Newton's coupling parameter G results in proportional evolutions for both 
the radial size and frequencies of the lunar motion. Slightly different orbital changes occur 
when a torque (indicated by L) acts on the orbit 

f 
r 

During the earlier years of the LLR experiment it has been the mean orbital radius signal 

Sr(t) me = - ^ r(t- t ) 

which has been used to measure G. But this involved estimating and subtracting a contri- 
bution to f which results from the orbital torque exerted on the Moon by the ocean tidal 
bulges on Earth which, because of friction, lag in angle from the direction toward the Moon. 
The inclination and 18.6 year precession of the lunar orbit's plane result in a modulation of 
the tidal contribution to f which helps the separation of the two perturbations after accu- 
mulation of sufficient years of data. But the data set produced by LLR has in recent years 
become sufficiently extended in time so that the range signals associated with frequency 
shifts, which grow quadratically in time, are becoming dominant in the fit for G. Recall that 
the two lunar phases can be expanded in terms of initial phase, rate, and acceleration 

D(t) = D + D(t-t ) + ±D(t-t ) 2 + ... (15) 



G L 
-G +2 ~L 

G L 



A(t) = A + A (t - Q + X -A (t - t f + 



(16) 
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The synodic frequency is by definition equal to difference of lunar sidereal rate and Sun's 
sidereal rate around the Earth 

D = cu - tt 

while the Moon's anomalistic rate is derivable from the underlying equation of motion and 
can be expressed in the form 

• 3tt 2 225 Q 3 , . .GM 

A = u - ... - 7 + l/2 — + ... 

4w 32 uj l c l R 

which consists of the classical Newtonian expression plus relativistic modifications, with the 
dominant geodetic precession contribution shown. From these two expressions the solar 
sidereal rate and its acceleration can then be expressed 

\ rS 3 (A-D) 2 

Q = A-D + -- : — — + ... 

4 A 

Q = A-D + ... 

While the lunar phases A and D suffer accelerations due to any tidal torques acting between 
Earth and Moon, the solar rate Q is not affected by the tidal torques. Acceleration of this 
solar rate is therefore a rather pure measure of a time variation of G. Noting from Equations 
(15-16) that the partials for A and D will grow in amplitude quadratic in time 

° R '"'" \t 2 (L ecc cos(A) — L evc cos(2D — A)) (17) 



OA 2 

dRcalc 



dD 



r (L var cos(2D) + L evc cos(2D — A)) (18) 



it follows that the formal error in measuring G decreases as the inverse square of the time 
span T of LLR observations. For a uniform time distribution of observations one obtains 



G in (5G\ /360 1 a 

with — 



G 2fi \GJ RMS ]j N QT 2 M L 2 + 3L 2 



with a being the rms size of individual range measurement errors, and N is the total number 
of measurements spread over the time T. A recent fit of almost 30 years of LLR data yields 
the excellent measurement constraint [10] 

| = (0 ± 1.1) 10- 12 (19) 

This amounts to about 1/60 the observed Hubble expansion rate of the universe. And with 
this measurement precision now growing quadratic in time, LLR should continue indefinately 
being at the cutting edge in supplying the measurement of G. 
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4. Additional Yukawa Interaction? 



When the supplementary interaction given by Eq. (7) is of a Yukawa nature, /i ^ 0, it 
contributes to precession of periastron for a near-circular orbit of radius r in amount 

S(cu — uj ) 1 KiKj 2 



2 GMiMj 



(jjr) exp(—/j,r) 



with w and uo Q being the orbit's sidereal and eccentric frequencies, respectively. And this 
perturbation of precession rate also occurs if the Yukawa force is metric, Ki ~ Mj, or non- 
metric. With the Moon's perigee precession rate measured to precision .07 mas/y and 
showing no anomaly, then for Yukawa ranges in the vicinity of that for maximum sensitivity 
of lunar perturbation, fir = 2, the strength of the Yukawa force is decisively constrained 

\K e K m \ ^ 5 1Q _ 12 / 4 _ 2) \ 



GM e M m ~ 



5. Gravitomagnetism. 



Line C of the complete N-body gravitational equation of motion given by Eq. (1) indicates 
a post-Newtonian gravitational force proportional to the velocities of both bodies in the in- 
teraction, and in analogy with electromagnetic theory, it has been called the gravitomagnetic 
interaction. From line C of Eq. (1), this acceleration is 

dux ■ 

5Si = (2 + 27) TT fa Vi ■ vj - fij ■ Vi Vj) 

It often has been claimed that the presence of gravitomagnetism within the total gravita- 
tional interaction has not been experimentally confirmed and measured. Indeed, different 
experiments have been under development to explicitly observe the effects of this historically 
interesting prediction of general relativity. But this gravitomagnetic acceleration already 
plays a large role in producing the final shape of the lunar orbit, albeit in conjunction with 
the rest of the total equation of motion; the precision fit of the LLR data indicates that 
gravitomagnetism's presence and specific strength in the equation of motion can hardly be 
in doubt. Because both the Earth and Moon are moving in the solar system barycentric 
frame — the frame in which the dynamical equations are formulated and then integrated 
into orbits — a gravitomagnetic interaction exists between these two bodies, the Earth hav- 
ing velocity V(t) and the Moon's being V(t) + u(t), as seen in Figure (2). As a result of 
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these mutual motions, perturbations to the Earth-Moon range from the gravitomagnetic 
acceleration result proportional to both V 2 and Vu 

Gm I 4 V 2 2 V11 \ 

- 530 cos(2L>) + 525 cos(L>) cm (20) 

As previously discussed, the amplitudes of the lunar motion at both these periods (monthly 
and semi-monthly) are determined to better than half a centimeter precision in the total 
orbital fit to the LLR data. It would be impossible to understand this fit of the LLR data 
without the participation of the gravitomagnetic interaction in the underlying model, and 
with strength very close to that provided by general relativity, 7 = 1. As in electromagnetic 
theory, the velocity-dependent force terms on lines C and D of Eq. (6) can individually be 
changed by formulating the dynamics in different frames of reference, but the very ability 
to reformulate the equations of motion in different frames without introducing new frame- 
dependent terms depends on the local Lorentz invariance (LLI) of gravity. It is the entire 
package of velocity-dependent, post-Newtonian terms which includes the gravitomagnetic 
terms, lines C plus D of Eq. (6), that produces the LLI; the Eddington parameter 7 represents 
the only freedom in the structure of this LLI package. Our confidence in the exhibited 
structure of this total collection of velocity-dependent terms is established in proportion to 
the precision with which the various preferred frame, LLI- violating effects in the solar system 
proportional to W 2 , WV, and Wu have been found to be absent [5]. LLR has been one 
of the main contributors in establishing gravity's LLI through null measurements of several 
VF-dependent effects [6, 4, 18, 7]. 



6. Inductive Inertial Forces. 



Inductive forces are shown on line E of Eq. (1); in such forces the acceleration of one mass 
element induces an acceleration of another proximite mass element (e.g., % and j in Figure 
2). From line E of Eq. (1) we have 

SS t = £ ((4 7 + 3)a, + a, • f y r y ) (21) 

These accelerations play a key part in altering the inertial masses of the Earth and Moon 
because of their internal gravitational binding energies; either the absence or an anomalous 
strength of these inductive forces would translate directly into differences between the ac- 
celeration rates of these whole bodies toward the Sun. A polarization of the Moon's orbit 
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in the solar direction, as previously discussed, would result. The forces, Eq. (21), acting 
between the mass elements of Earth, for example, by themselves would lead to an anomalous 
polarization of the lunar orbit of very large magnitude 



5r(t); = 130 cos(£>) meters 



(22) 



Only when these inductive forces are combined with the other post-Newtonian inertial forces 
shown on line F of Eq. (1) does the total inertial self force of a body become 



Sf = 




G . mirrij 



ho 



rij 



i 



G x mirrij 



1,3 



r 3 r *i 

ij 



The first line of this total self force is now the expected inertial force due to the internal 
kinetic energy and gravitational binding energy within the body. The second line represents 
contributions to the body's internal virial which, when totaled over all internal force fields, 
vanishes for a body in internal equilibrium and experiencing negligible external tidal-like 
forces. These self forces of a body are an integral part of the determination of the total 
gravitational to inertial mass ratio of bodies discussed previously, and in general relativity 
they are cancelled by equal contributions of internal energies to a body's gravitational mass. 
They were explicitly discussed here in order to show the large size of such inductive force 
contributions which must necessarily be taken into account in the fit of theory to the LLR 
data. 
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